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ABSTRACT 

We use a robust sample of 11 z ^ 7 galaxies (zsso-dropouts) to estimate the stellar mass density 
of the universe when it was only ^750 Myr old. We combine the very deep optical to near-Infrared 
photometry from the HST ACS and NICMOS cameras with mid-Infrared Spitzer IRAC imaging avail- 
able through the GOODS program. After carefully removing the flux from contaminating foreground 
sources we have obtained reliable photometry in the 3.6 jiva and 4.5 ^m IRAC channels. The spectral 
shapes of these sources, including their rest frame optical c olors, strongly support the ir being at z^7 
with a mean photometric redshift of (z) =7.2±0.5. We use iBruzual fc CharlotI (j2003f ) synthetic stel- 
lar population models to constrain their stellar masses and star formation histories. We find stellar 
masses that range over 0.1 - 12x 10^ Mq and average ages from 20 Myr to up to 425 Myr with a mean 
of ~ 300 Myr, suggesting that in some of these galaxies most of the stars were formed at z > 8 (and 
probably at z > 10). The best fits to the observed SEDs are consistent with little or no dust extinction, 
in agreement with recent results at z ~ 4 — 8. The star formation rates (SFR) are in the range from 
5-20 Mq yr~^. From this sample we measure a stellar mass density of 6.6^3 3 x 10^ Mq Mpc~'^ to a 
limit of Muv,AB < —20 (or 0.4L*^3). Combined with a fiducial lower limit for their ages (80 Myr) 
this implies a maximum SFR density of 0.008 Mq yr~^Mpc~'^. This is well below the critical level 
needed to reionizc the universe at z ~ 8 using standard assumptions. However, this result is based 
on luminous sources (>L*) and does not include the dominant contribution of the fainter galaxies. 
Strikingly, we find that the specific SFR is constant from z ~ 7 to z ~ 2 but drops substantially at 
more recent times. 

Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

Direct observations of galaxies at high redshift 
from complete, well-defined searches place strong con- 
straints on galaxy formation and evolution mod- 
els. Extensive studies have been made of galax- 
ies selected by th e dropout technique out to red- 
shift 6 (see e.g.. iStanwav et al] 120031: iBunker et al.l 



200llYan & W indhorst ll2004tlBouwens et al.ll200(ill2007r 
McLure et al I2OO9 a) thanks in part to the capabilities of 
the HST ACS and NICMOS cameras. Studies at higher 
redshifts, however, have been much more challenging. In 
particular, one of the key issues has been the derivation 
of stellar masses. These masses can provide both strong 
additional constraints on formation models and informa- 
tion about the star formation rates at even earlier times, 
effectively opening a window towards the earliest phases 
of galaxy formation. Deriving masses at higher redshifts 
does however face some distinct challenges. 

At redshifts z > 4, ACS and NICMOS only access 
the rest-frame UV continuum, which is not a reliable 
tracer of stellar mass. To obtain better constraints on 
this quantity through SED fits it is necessary to ex- 
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tend the observations to the rest-frame optical. Although 
Spitzer data presents challenges due to its large instru- 
mental PSF, these rest-frame optical measurements can 
be made at high redshifts with the 3.6 tiui and 4 .5 /xm 
IR channels of the IRAC camera (jFazio et al.ll2004f ). The 
IRAC 3.6 ^m and 4.5 ^m bands probe the rest-frame op- 
tical fluxes around 0.5 — 0.6 ^m of z ~ 5-|- sources and 
reach to >26 AB mag at 5a in deep >23 hr integrations. 
This is a remarkable achievement for a 0.8-m telescope. 
Spitzer will continue to play a unique role in the deter- 
mination of fundamental properties like ages and stellar 
masses of the earliest galaxies until JWST is launched. 
These Spitzer and HST data have permitted estimates 
of stella r masses for large numbers (>150) of z ~5-6 
sources (lYan et all 120051: lEvles et al.ll2 005>: 'Egam i et al" 
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Yan et al.ll2006t iStark et al.ll2009 : Mobashe r et al 



Bradlev et al.l 120081 : iZheng et al.l 120091) . One sur 



prising early finding was the number of qui te massive 
^ 10^° Mq galaxies in z ~ 6 samples fe.g lYan et al.l 
[200l [E^desetal] 12005?) . It suggested that z ~ 5 - 6 
galaxies exhibited substantial amounts of star formation 
at much higher redshif ts and ear lier times, well into the 
epoch of reionization (Stark et al. 2007; Yan et al. 2006). 

These analyses have also opened the possibility of 
looking at the SFR-Mass relation at higher redshifts, 
or equivalently, at the specific SFR or SSFR (i.e., 
SFR/Mass). The SSFR tells us how fast galaxies are 
growing with respect to their current stellar masses. 
While at low redshifts (0 < z < 1) it has been shown 
that g a.laxies grew faster in earlier times than they do 
today (jNoeske et al.ll200"7l ). it has been suggested that 
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at higher redshift (4 < z < 6 ) the SFR-Mass relation 
remains much more constant (iStark et alj 120091) . This 
transition in the SSFR is an interesting and important 
result and needs to be assessed over a long redshift base- 
line. 

It is imperative then to try to extend these studies 
to z > 7 and investigate what are the characteristics of 
the probable progenitors of these rather massive z ^ 6 
sources. Doing so, however, has been challenging due 
to the difficulty in obtaining deep enough near-IR and 
optical data to robustly identify z ~ 7 sources. Such 
sources only start to become "common" at >26.5 AB 
mag. Thus, early selections of z ~ 7 z-dropouts from 
the HUDF only included a handful of c andidates, and 
the situatio n has only i mproved slowly (jBouwens et al.l 
[2004c, 200l lOesch etal. 2009^. Consequently, very lit- 
tle has been publ ished on the stellar masses of z ~ 7 z - 
dropout galaxies (lEgami et al.ll2005t ILabbe et"aI1l2006D . 
iLabbe et al.l (|2006[ ) performed stellar population model- 
ing of four ACS zgso-band dropout sources found in the 
HUDF, and derived stellar masses of - 10^ - 10^° Mq, 
ages of ^--^100-200 Myr, and was able to estimate a stellar 
mass density of 1.6 x 10^ Mq Mpc^^ (to 0.3L*^3). 

Fortunately, as a result of continued efforts to se- 
lect z ~ 7 galaxies from the growing quantity of deep 
near-IR NICMOS data, Bouwcns ct al. (2009, in prepa- 
ration) have succeeded in substantially expanding the 
size of current z ~ 7 z-dropout selections, and now 
14 NICMOS-sclcctcd, rather luminous z ~ 7 z-dropout 
galaxies are known that are amenable to stellar mass 
estimates. The recent advent of WFC3 on HST has 
expanded the number of known z ~ 7 sources and 
now ~25 sources are known at z 7 (lOesch et al.l 



2p09a'; 'Bouwens et all l2009bl: iMcLure et al. 

Bunker et^ 20091: iCastellano et al.l 120091: lYan et al 



2009bl: 



20091: lOuchi et al.ll2009l : IHicke7 et al 1l2009HWilki~ etal 



20091) . Most of the newly discovered galaxies, how- 
ever, are too faint to attempt stellar mass estimations 
on individual basis, although valuable information can 
be obt ained from these sam ples through stacking analy- 
sis (see ILabbe et al.l 120 Kfi ). Galaxies in the Bouwens ct 
al. (2009, in preparation) sample extend from ^25.5 AB 
mag to ~27.8 AB mag and are found in the HUDF, and 
in and around the two wide-area GOODS fields. These 
candidates also possess very deep >23.3 hr coverage with 
IRAC from the Spitzer GOODS program and present us 
with a unique opportunity to better understand what 
the typical properties (ages, stellar masses) of z ~ 7 z- 
dropouts are. 

Here we take advantage of the larger sample of 14 zgsp- 
dropout sources identified by Bouwens et al. (2009, in 
preparation) to estimate the typical properties of z ^ 7 
galaxies. The much larger size of current samples al- 
lows us to get a better handle on the typical properties 
of z ^ 7 galaxies than was possible from the smaller 
and brighter samples previously available (jLabbe et al.l 
|2006[) . For example, sizeable variations in the M/L ra- 
tios of individual galaxies can considerably skew the av- 
erages for the population as a whole. The present sample 



^ The ILabbe et al.l 1)20101 ) work focuses exclusively on the ultra- 
faint, sub-L* sources. Only sources UDF-387-1125 and UDF-3244- 
4727 from this work (our faintest z gsp— dropouts in the UDF) are 
also included in ILabbe et al.l 1)20101 ). 



also include galaxies from two independent lines of sight, 
i.e., the HDF-North GOODS and CDF-South GOODS, 
so the results should be much less impacted by cosmic 
variance and more representative of the cosmic average. 
In estimating the rest-frame optical fluxes for our z ^ 7 
candidates from the IRAC, we will take adv antage of the 
well-te sted deble nding process descri bed by ILabbe et al.l 
()2006[) (see also IWuvts et abl I2007D which enables us 
to estimate fluxes when there is moderate overlap with 
nearby sources. This also allows for a larger sample, since 
the well-known problems with blending and confusion in 
IRAC data have limited previous studies. 

We provide a brief outline for this paper here. We 
present the sample selection and observational data we 
have used in this work in §2 and §3 respectively. In 
§4, we will describe the photometry of the sample with 
particular emphasis in the deblendcd photometry from 
the IRAC channels. We devote §5 to a discussion of 
the photometric redshifts obtained and possible contam- 
ination and in §6 we present the procedure and results 
in the process of fltting synthetic stellar populations to 
the observed SEDs (including our SSFR results). §7 
presents our estimated stellar mass density (SMD) at 
z ^ 7 and confldence intervals. We compare our re- 
sults to previous work in §8 and present a discussion of 
them in §9. We summarize our conclusions in §10. All 
magnitudes quoted in the paper are in the AB system 
(jOke fc Gunnl[T983l ) . We have used cosmological param- 
eters = 0.7,17a/ = 0.3 and Hq = 70 kms s^^ Mpc^^ 
to facilitate comparison to previous works. 

2. SAMPLE SELECTION 

The present selection of z ~ 7 z-dropout candidates is 
based upon ^ 80 arcmin^ of very deep optical ACS, near- 
IR NICMOS, and IRAC data available over and around 
the two GOODS fields - including the HUDF and HDF- 
North (see §3 for a more detailed description of these 
data). This selection is described in Bouwens ct al. (2009, 
in preparation). 

Candidates z ^ 7 galaxies were required to satisfy a 
two-color Lyman-Break Galaxy (LBG) criterion adapted 
to z ^ 7 - showing a strong z — J break but possessing 
a blue J — H color (redward of the break). In detail, the 
color criterion used is: 



A(; 



■850 



(^850 — ^iio)ab > 0.8 
JllOlAB > 0.8 + 0.4( Jiio - Hi6o)ab 



where A represents the logical AND operation and where 
the colors were measured in a '-^0.4"-diameter aper- 
ture using SExtractor in double image mode with the 
iJieo band image being the detection image. Candidate 
z > 7 sources were also required to be undetected to less 
than 2a in the optical ACS bands (i3435, Vijoe, and 1775). 
Sources were also removed from the sample if they were 
detected at l.Str in two or more of the optical bands. 
Only detections of at least 5(t in the iJieo-band images 
(measured in apertures of 0.6" in diameter) were consid- 
ered to ensure that our candidates corresponded to real 
sources. The candidates -with GNS IDs- were also re- 
quired to have have (iJieo — 5.8 /xm) colors bluer than 
2.5 (5.8 /xm photometry measured in 2"diameter aper- 
tures). This provided added confidence that the sources 
were not low-redshift interlopers. 
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TABLE 1 
Summary of available data. 



Filter 



Depth Reference 
[la] 



HUDF 



ACS - B435 


29.9 


1 




ACS - i;606 


30.4 


1 




ACS - 2775 


29.9 


1 




ACS - 2850 


29.4 


1 




NICMOS - Jiio 


28.7 


5, 6] 


NICMOS - HiQo 


28.3 


5, 6 


ISAAC - Ks 


27.3 


11 




IRAC - 3.6/im 


27.5 


15 




IRAC - 4.5/im 


26.8 


15 




IRAC - 5.8/im 


25.8 


15 




IRAC - 8.0/im 


25.7 


15 





GOODS South 



ACS - B435 


28.8 


2, 


3 




ACS - vaoa 


29.1 


2, 


3 




ACS - i775 


28.5 


2, 


3 




ACS - ^850 


28.3 


2, 


3 




NICMOS - Jiio 


28.0 


7, 


8 




NICMOS - Hie.0 


27.9 


7, 8 




IRAC - 3.6/im 


27.2 


15 






IRAC - 4.5/im 


26.4 


15 






IRAC - 5.8/im 


25.4 


15 






IRAC - 9,.0^lm 


25.3 


15 






GOODS North 








ACS - S435 


28.7 


2, 


3 




ACS - veoe 


29.0 


2, 


3 




ACS - i775 


28.6 


2, 


3 




ACS - 2850 


28.5 


2, 


3 




NICMOS - Jiio 


27.8 


7, 


8 




NICMOS - Hieo 


27.7 


7, 


8 




MOIRCS - Js 


27.1'' 


12, 13, 14] 


MOIRCS - Ks 


27.4=^ 


12, 13, 14 


IRAC - 3.6/im 


27.2 


15 






IRAC - A.bfim 


26.4 


15 






IRAC - 5.8/im 


25.4 


15 






IRAC - 8.0/im 


25.3 


15 






HDFN'" 








ACS - S435 


29.2 


2, 3] 


WFPC2 - S^^go 


30.2 


41 






ACS - 1)606 


29.4 


2, 3] 


WFPC2 - 

ACS - 2775 


29.0 


41 






28.6 


2, 31 


WFPC2 - /814 


28.8 


41 






ACS - Z850 


28.23 


2, 31 


NICMOS - Jiio 


27.9 


9, 


10] 


NICMOS - Hiao 


28.1 


9, 


10] 


MOIRCS - Ks 


27.6 


12, 13, 141 


IRAC - 3.6/im 


26.8 


15 






IRAC - A.bfim 


26.2 


15 






IRAC - 5.8/im 


25.8 


15 






IRAC - 8.0/im 


25.7 


15 







Re feren ces. — (1) [B cckwith et al. 
200^: (2) IGiavalisco ct al. 2004; (3) 
Bouwens et al. 2007; (4) Williams ct al.l 
19961: (5) , Thonipson ct al, 2005, ;; (6) 
Oesch et al.l I2009bl i (7) Conselice et 
al. 2009, in preparation; (8) Bouwens 

et al. 2009, in prepara tion; (9) 

Thompson et al.l [19991: ( 10) IDickinsonI 
19991: (111 ILabbe et all 120031: (12) 
Kaiisawa etall 12003: (13) [Duchi et al.l 
2001 (14) IBouwens et al.l 120081 : (15) 
Dickinson et al. 2009, in preparation 



Note. — The depths of the data in the 
different fields. Optical to NIR were esti- 
mated using 0.9" diameter apertures and 
corrected to total assuming stellar pro- 
files (aperture corrections of 5% and 20% 
for the ACS and NICMOS data respec- 
tively). Apertures of 1.8" were dropped 
randomly in empty regions of the IRAC 
images and then the fluxes corrected to 
total assuming PSF profiles (aperture 
correction multiplicative factors of 2.4, 
2.7, 3.5, and 3.7 for the 3.6 iim , 4.5 
um , 5.8 um , and 8.0 um images re- 



By applying the above criteria to ^80 arcmin^ of deep 
NICMOS data, Bouwens et al. (2009, in preparation) 
identified 14 ^ ~ 7 ^-dropout candidates (sec Table [2]). 
The candidates range in iJi6o,AB magnitude from 25.5 
to 27.6 AB mag and are typically only marginally re- 
solved at NICM0S/NIC3 resolution (FWHM of the PSF 
is ^ 0.37"). Contamination rates for the sample are de- 
termined by running a number of photometric scattering 
simulations on a fake sample of low redshift sources. This 
sample is constructed to match the color distribution of 
observed galaxies in the 24.5 <i/i6o< 26 range. We have 
added noise to the fluxes in each individual band accord- 
ing to the depths of the different fields. By applying 
the selection criteria previously described, we character- 
ize when and with what frequency low-redshift contam- 
inants enter our selection. A more detailed discussion 
of the procedure is given in IBouwens et"all (|2008D and 
Bouwens et al. (2009, in preparation). Through this tech- 
nique we find that the contamination for the sample (by 
lower redshift galaxies, time-variable sources, low-mass 
stars, for example) is expected to be just ~10% of the 
sample for sources found over the HUDF and HUDF05 
fields, and only ^20% for sources within the GOODS 
fields. 

We only use those 11 z-dropout candidates from the 
Bouwens et al. (2009, in preparation) selection that have 
deep (23-1- hr) IRAC data in both 3.6 and 4.5 /xm and 
which were in the Bouwens et al. (2009, in prep) z- 
dropout sample as of June 2009. 

3. OBSERVATIONAL DATA 

The very deep optical, near-IR, and mid-IR IRAC data 
available for our z ~7 2-dropout candidates permit us to 
study the properties of these sources in great detail. A 
summary of the available imaging data for the candidates 
we have in our various search fields is given in Table [1] 

The deep near-IR data we have available for our can- 
didates at ~^1.1 jiva and ~^1.6 ^m comes from NICMOS. 
For the candidates within the GOODS fields, the NIC- 
MOS data reach to depths of --26.8 AB mag and --26.7 
AB mag (5(t, aperture flux within a 0.6"— diameter, used 
for detection ) in the Jno and i?i6o bands, respectively. 
For purposes of the SED fitting we use a larger (0.9"- 
diameter) aperture to minimize differences in the aper- 
ture corrections among the optical, NIR and mid-IR im- 
ages. The Icr limits relevant for the SED fitting are 
~ 27.8 and ^ 27.9 AB mag (Jno and i/ieo bands, cor- 
rected to total fluxes using aperture corrections of ^ 20% 
derived from stellar profiles). The near-IR data over the 
HUDF reaches some ~0.5 AB mag deeper. The FWHM 
of the NICMOS PSF is 0.34" and 0.37" in the Jno and 
HiQQ bands, respectively. 

We also have very deep ~2.2 ^m K-band data avail- 
able for our 2; ~ 7 z-dropout candidates over the HUDF 
and in the central region of the HDF-North. These data 
are particularly valuable for providing a constraint on 
the UV-continuum slopes. The i^-band data over the 
HUDF correspond to the 40 hr integration in the best 
seeing conditions in the K^-band filter at ISAAC (VLT) 
and PANIC (Magellan) (jLabbe et al.ll2006D . We estimate 
Icr depths of ~ 27.1 - 27.4 AB magnitudes in 0.9" di- 
ameter apertures. The other two sources are in the 
HDFN and their K-band data come fro m the deep Sub- 
aru GTO MOIRCS imaging campaign (|Kaiisawa et al.1 
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TABLE 2 
Summary of Photometry. 



ID 


Field 


^43 5 




*775 


2850 


-'no 






3.6/im 


4.5/im 


UDF-640-1417 


HUDF 


>29.9 


>30.4 


>29.9 


27.8±0.2 


26.5±0.2 


26.1±0.2 


25.8±0.2 


25.2±0.1 


25.9±0.4 


UDF-983-964 


HUDF 


>29.9 


>30.4 


>29.8 


>29.2 


26.7±0.2 


26.8±0.3 


26.6±0.6 


26.2±0.3 


26.6±0.9 


UDF-387-1125 


HUDF 


>30.0 


>30.6 


>30.0 


29.1±0.8 


27.7±0.4 


26.9±0.3 


>27.4 


>27.5 


>26.8 


UDF-3244-4727 


HUDF 


>30.0 


>30.4 


>30.0 


>29.4 


27.6±0.3 


26.9±0.3 


>27.1 


26.5±0.4 


26.7±1.0 


GNS-zDl 


GOODS-S 


>28.6 


>29.0 


>28.3 


>27.9 


26.5±0.3 


26.2±0.2 




24.9±0.1 


25.4±0.4 


GNS-ZD2 


GOODS-S 


>28.9 


29.0±1.0 


>28.4 


>28.4 


27.0±0.4 


26.7±0.3 




26.4±0.6* 


26.2±0.9 


GNS-zDS 


GOODS-S 


>29.1 


>29.4 


>29.0 


>28.7 


27.3±0.7 


26.7±0.3 




25.7±0.3 


25.9±0.7 


GNS-ZD4 


GOODS-N 


>28.5 


>28.8 


>28.4 


>28.4 


26.7±0.4 


26.3±0.3 




25.2±0.2 


25.4±0.4 


GNS-zDS'' 


GOODS-N 


>28.8 


>29.1 


>28.7 


28.1±0.7 


26.2±0.3 


25.3±0.1 


25.4±0.2 


24.6±0.1* 


25.4±0.3 


CDFS-3225-4627 


GOODS-S 


>28.7 


>29.0 


>28.6 


>28.4 


26.8±0.4 


26.4±0.3 




26.1±0.3 


26.0±0.5 


HDFN-3654-1216'^ 


HDFN 


>29.2 


>29.4 


>28.6 


26.5±0.2 


26.2±0.2 


25.9±0.2 


26.4±0.4 


25.2±0.3* 


25.7±0.7 


Mean SED^ 




>30.3 


>30.6 


>30.1 


28.41±0.29 


26.70±0.09 


26.27±0.07 


26.68±0.30 


25.58±0.07 


25.91±0.18 



Note. — Magnitudes are total and in the AB system. Optical to near infrared photometry measured in 0.9" diameter apertures with aperture 
corrections of 5% - 20% (derived assuming stellar profiles). Mid infrared IRAC photometry performed on "cleaned" images with 1.8" aperture 
diameters. Aperture correction factors in this case are 2.4 and 2.7 in the 3.6 and 4.5 fim channels respectively. Upper limits and error bars are 1 
fj. 

^ K-band from either MOIRCS or ISAAC depending on the Field, see Table [T] 

^ MOIRCS is band imaging is deep enough at the location of this source. We measure Js=25.6±0.3. 
There are additional WFPC2 optical imaging constraints for HDFN-3654-1216: B> 30.2, V> 29.0, i >28.8. 

This mean SED was constructed rescaling all the SEDs so that the _ffi60-band fluxes coincide with the mean value and then averaging all the 
other bands. The 3.6fj,m IRAC fluxes from sources GNS-zD2, GNS-zDS, and HDFN-3654-1216 were excluded from the average. See section §6.2 
and figure |4] 

After subtracting a bright nearby neighboring from the S.Gfim image of this source, sizable residuals remain, and so the quoted 3.6fim fiux 
measurements for these sources may suffer from systematic errors. 

[2006t lOuchi et al.l [2007t see iBouwens et all [20081 for a 

description of our reductions). We have only made use 
of the deepest GTO pointing which reaches down to 25.4 
total AB magnitudes at Scr (0.9" diameter apertures). 
The FWHM for the /-iT-band PSF is -0.5". 

Deep optical ^435, vqqg, *775: and zsso observations 
are available for our candidates with ACS and typically 
reach to depths of — 30 mag in the HUDF and ~1.5 
mag shallower in the rest of our fields. For the single z- 
dropout candidate in the HDF-North WFPC2 field, we 
have the very deep (> 28 mag at 5a) i?45o,V606, and 
/814 WFPC2 observations - which permit us to set very 
strong constraints on the strength of the Lyman Break. 

The Spitzcr IRAC imaging data from the GOODS pro- 
gram (Dickinson et al. 2009, in preparation) provide us 
with deep rest-frame optical coverage on our z ^ 7 z- 
dropout candidates - which is critical for estimates of 
the stellar mass in these sources. Two exposures of ^23- 
hr each were taken in two different epochs with the IRAC 
camera - rotated by 180 deg - and overlapping in the cen- 
ter of the GOODS field. The region of overlap contains 
the HUDF in the GOODS South. We find AB magnitude 
detection limits for point sources of 27.4, 26.6, 25.4, and 
25.3 for the 3.6, 4.5, 5.8, and 8.0 fim channels respec- 
tively (1 cr, measured on apertures of 1.8" in diameter 
and corrected to total flux assuming stellar profiles with 
aperture corrections of 2.4, 2.7, 3.5, 3.7 respectively - 
multiplicative factors). These limiting depths (for single 
epoch ^23.3 hr IRAC observations) were estimated by 
dropping apertures at random empty regions of the sky 
and measuring the flux variations. A good summary o f 
the IRAC observations is provided in lLabbe et al.l (j2006l ) 
and iStark et al.l (I2007D . for example. In this work we 
make use of the reductions of Data Release (DR) 3 of 
epoch 1 observations and DR2 of epoch 2 of the GOODS- 



S field. In the case of the GOODS-N field wc make use 
of the reductions of Data Release (DR) 2 of both epochs. 

4. PHOTOMETRY 

Optical/near- IR Photometry: Optical to NIR fluxes 
were measured in standard 0.9"-diameter circular aper- 
tures. We corrected these measured fluxes for the miss- 
ing light outside these apertures assuming stellar profiles. 
These latter corrections increased the measured fiux by 
5% - 20% depending on the band. All the sources in our 
sample were sufficiently separated from their neighbors 
that wc could use this simple approach, except for the one 
in the HDFN field where there was a faint but very close 
(almost overlapping) neighbor. This made it impossible 
to use the standard circular aperture to measure the flux 
of HDFN-3654-1216 (ID from Bouwens et al. 2009, in 
preparation) . To ensure our measurement of its flux was 
not contaminated by the flux of this neighbor , we fitted a 
PSF profile to the neighbor using GALFIT (|Peng et alJ 
|2002| ) and subtracted it from the image before measuring 
the flux of this z-dropout. 

IRAC Photometry: One of the biggest challenges in es- 
timating the masses of our z ^ 1 candidates is the acqui- 
sition of reliable mid-IR fluxes for these candidates from 
the available IRAC data. The extremely broad PSF of 
this instrument and large pixel sizes make the images ex- 
tremely crowded and so fluxes from neighboring sources 
spill over onto each other. To overcome this issue, a wide 
variety of different approaches have been developed, al- 
most all of which involve modeling the IRAC image with 
a number of smoothed sources of varying flux. In the 
most common cases, the model light profiles are theoret- 
ical like the ones produced by GALFIT. The use of these 
models can result in systematic errors in the photometry 
if the sources have irregular or clumpy spatial profiles. 
We have used here the technique described in 
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iLabbe et al.l (|2006( ). which consists in the creation of 
an empirical hght profile based in the higher resolution 
NICMOS images. This tec hnique uses a segrnentat ion 
map created by SExtractor (jBertin fc Arnoutill996f ) to 
define the boundaries of each source in the area to be 
cleaned (wc use a 2cr threshold to ensure all the possibly 
relevant neighbors are fitted) and use the light profiles 
of the sources within those boundaries as the empirical 
light profiles (assuming a similarity between the profiles 
at both 1.6 fim and 3.6 fim). The individual profiles are 
then convolved with a carefully constructed kernel (based 
on the instrumental PSFs) to simulate how they would 
look like in the IRAC images (modulo a normalization 
factor). Finally we fit for the total flux of each neighbor 
and subtract them off the image. Instead of using the 
z-dropout flux measurement determined from these fits, 
wc subtract off the flux from the neighbors and perform 
standard aperture photometry in relatively small aper- 
tures on the "cleaned" image. We find that the optimal 
aperture diameter for maximizing the S/N of our flux 
measurements (assuming point sources) is ~1.8" in all 
the channels. The aperture corrections are 2.4, 2.7, 3.5, 
and 3.7 (multiplicative factors) in the 3.6, 4.5, 5.8, and 
8.0 /im channels respectively. The errors in the measured 
flux include both the typical variations on the sky bright- 
ness, and the uncertainty in the flux removed from the 
aperture. 

A simple inspection of our images shows that 5 of the 
11 z ^ 7 z-dropout candidates in our sample (i.e., GNS- 
zD2, GNS-ZD5, UDF-387-1125, UDF-3244-4727, HDFN- 
3654-1216) are severely blended with bright foreground 
sources. The other six sources are somewhat more iso- 
lated but with the size of the instrumental PSF it is ob- 
vious that flux from neighboring sources will spill over 
onto these candidates (often contributing >20% of the 
light within the 1.8"-diameter aperture centered on our 
candidates) . From the 5 severely blended sources we are 
able to satisfactorily "clean" the images in at least 2 cases 
(we will discuss the other three cases GNS-zD2, GNS-zD5 
and HDFN-3654-1216 in the next paragraph). Wc could 
check the consistency of the method for 5 of the sources 
that were imaged in both epochs. As was mentioned 
before, there is a 180 deg rotation in the IRAC camera 
(and thus in the asymmetric instrumental PSF) between 
the two epochs, which makes the light profile models (for 
the neighboring sources) almost independent. Obtaining 
consistent fluxes in these two images is a good indication 
of the reliability of the method. For the sources where 
IRAC observations are from >1 epoch (and consistent 
within the measurement errors), we average the fluxes 
and combine the errors accordingly which reduces them 
by a factor \/2. 

In the cases of GNS-zD2 and GNS-zD5 strong color 
gradients in the closest neighbor cause the model light 
profiles to be inadequate and so large residuals are ev- 
ident (after subtracting the flux from the neighbors). 
GNS-zD5 was imaged in both epochs and we find a dif- 
ference of 2(7 between the measurements in the 3.6 fj.m 
channel but a much cleaner residual (and better agree- 
ment) in the 4.5 /xm image. A similar discrepancy was 
found for the 3.6 fim flux measurement for HDFN-3654- 
1216. In this case the poor subtraction is not due to 
strong color gradients but to the extreme proximity of 



its neighbor. In all three cases we adopt the single epoch 
uncertainty. Given the small size of our present sample 
of z-dropouts, we chose to keep these three sources in the 
sample - keeping in mind the caveat that the flux mea- 
surements for these sources could possess large system- 
atic errors. In summary, we are able to perform reliable 
cleaned photometry on 8 out of the 11 sources (73%). 

We obtain > 2a detections for 9 of the 11 sources in 
the 3.6 fj,m image and for 5 sources in the 4.5 fj,m im- 
age. The two sources with quite marginal (< 2a) detec- 
tions in both IRAC images are UDF-387-1125 and GNS- 
zD2. A simple stacking of the images of these two sources 
(adding both sources in both IRAC images) shows a sig- 
nificant detection, which provides evidence for the reality 
of the sources. 

Optical to mid-Infrared image stamps 4" x 4") for 
all the sources are presented in figure [TJ The 5.8 yum 
and 8.0 fj.m channels have been omitted because none of 
the sources are detected in those bands (as expected). 
The two epochs have been coadded when available. The 
measured magnitudes are summarized in table [2j 

5. PHOTOMETRIC REDSHIFTS 

Perhaps the most fundamental quantity that we can 
estimate for galaxies in our zgso-dropout sample is their 
redshift. We first explore the probable rcdshifts of 
our sources usiiig the photometric redshift code EAZY 
(jBrammer et al.ll2008[) . The code works by comparing 
the observed photometry with that predicted on the ba- 
sis of the specific SED templates. We use the default 
template set, which wa s derived from the Pegase popula- 
tion synthesis models (jFioc &: Rocca-Volmerangj|1997[) 
and optimized to reproduce the properties of galaxies in 
the range < z < 4. Comparisons of the photomet- 
ric redshift estimates output from EAZY with some of 
the deepest spectroscopic surveys available show minimal 
systematic errors and a scatter in Zspec — Zvhot of (7=0.034 
over the range < z < 4 (see table 1 in iBrammer et ahl 
l2008f ). The few sources at higher rcdshifts (z '--^ 6) with 
available spectroscopy are also in good agreement with 
its photometric redshift estimates. One relevant advan- 
tage over other photometric redshift codes is that EAZY 
works with fluxes instead of magnitudes and naturally 
handles negative measurements which are common in our 
optical bands. 

We use the available flux measurements in the i?435; 
Veoe, i775, ^850, -^iio, -f^ieo, 3.6 ^m, 4.5 ^m bands, and 
also the -B450, Vgoe, Hui K bands if available for our 
photometric redshift estimates. We do not include the 
IRAC 5.8 /im and 8.0 /im flux measurements of the 
sources (consistent with no detection) in the compari- 
son - since they do not help us to meaningfully discrim- 
inate between the competing redshift solutions. We re- 
strict the redshift range of our fits to z ^ 4—11 and 
adopt no redshift prior. Solutions at z ~ 1.5 are also 
possible for most of the sources but at lower probability 
for all of them. However, quantifying this is challenging 
given our poor knowledge of the demographics of galax- 
ies at both z ~ 1.5 and z ~ 7. Furthermore, the exis- 
tence of synthetic solutions at any of these rcdshifts does 
not necessarily imply the existence of real galaxies with 
the observed properties. We have provided an indepen- 
dent (short) description of the estimate of the fraction of 
low-redshift contaminants in §2. A more detailed discus- 
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Fig. 1. — Postage stamps of the zgso— dropouts. Each stamp is 4"x4" in size (~ 21 kpc at z = 7). The dashed circles indicate the 
apertures used for photometry (0.9" for the optical to NIR and 1.8" in the IRAC bands). All the sources are undetected in £435, U6O61 
and 1775 and only UDF-640-1417 and HDFN-3654-1216 arc weakly detected in zgsQ. The latest two sources are estimated to be the lowest 
redshifts of the sample. 3.6 fim and 4.5 /im stamps show the sources after the flux from the nearby neighbors have been fit and removed 
(i.e., they are "cleaned" images), with the two epochs of IRAC data added together when available. The cleaning process in the case of 
sources GNS-zD5 and GNS-zD2 has left nearby residuals (seen in white) attributed to the close and bright nearby sources visible in the 
other bands. To a lesser degree this seems to also be the case for sources UDF-640-1417 and UDF-3244-4727. Since these latter sources 
were observed in both epochs, it was possible to check that the measured fluxes are consistent among them. As expected, all sources are 
undetected in the 5.8 and 8.0 fim channels. The sources have been placed in order of increasing Hiqq magnitude from top to bottom. 
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Fig. 2. — The photometric redshifts of our 2 ~ 7 z-dropout sam- 
ple (histograml . The soUd hne shows the redshift distribution that 
IBouwens et al.l lj2008l ) predict for the 2-dropout selection. The pho- 
tometric redshifts were estimated with EAZY (see §4). Typical un- 
certainties in the redshift for individual sources is Az ~ 0.6. The 
lowest redshift in the sample corresponds to object HDFN-3654- 
1216 at 2 = 6.2. This object presents the bluest zgso— Jiio color 
due to it being weakly detected in the zgso band image. 

sion is provided in the Bouwens et al. (2009, in prepara- 
tion) study where th e sample is described (but see also 
IBouwens et al.l [20081 which describes the estimate for a 
similar sample). The estimated contamination fraction 
for this sample is ~ 10% and ^ 20% for sources in the 
HUDF and outside the HUDF, respectively. 

We explore the solutions in the range 4 < z < 11 in 
a grid with steps of 0.01. At each step there is a x^{^) 
value associated with the best solution that is used to 
create a probability function p{z) cx cxp[— x^(z)/2]. The 
redshift is estimated marginalizing over this probability 
and probability contours are used to determine the 68% 
confidence intervals. Because of the shape of this func- 
tion at redshifts z > 7, the estimated redshifts are usually 
located somewhat above the absolute minimum x^- The 
typical uncertainties we estimate for individual sources 
are in the order of 0.6. Finally, the redshift distribution 
obtained from EAZY is z = 7.2 ± 0.5, similar to what 
would be expected based on the shape of the ACS zs5o 
passband and the color cuts imposed in the sample(sec 
figure [21). 

6. STELLAR POPULATION MODELING 

We can make sense of the present photometric selec- 
tion of z-dropouts in terms of their intrinsic properties 
- like age, stellar mass, or dust extinction - by mod- 
eling the spectral energy distributions observed for in- 
dividual sourc es. Such modeling i s now ubiquitous in 
the literature (iSawicki &: Yeel[T99l iBriiichmann &: Ellii 
[2000t [Parovich et al.l 120011: iLabbe et al.ll2007[ ). and has 
proven quite powerful in the interpretation of distant 
galaxies - both given the abundance (and quality) of 
photometric data and the plausibility of the photometric 
estimates. Indeed, studies have found reasonable agree- 
ment between stellar masses determined from such mod- 
eling: and those determined from the dynamics (|Erb et al.l 
[2009) . 

Here we model the stellar populations with the 



iBruzual fc CharlotI (|2003[ ) (BC03) spectral synthesis li- 
braries. Over the wavelength range covered by the 
present study, these libraries show almost no difference 
with respect to the newest Chariot fc Bruzual ( 2007) 
(CB07, e.g. IStark et al.l [200l ILabbe et al.l [20Tol) and 
have been shown t o be in reasona ble agreem ent with 
other libraries (e.g.. lMaraston|[2005[ ). We use a lSalpete"il 
(|1955[ ) Initial Mass Function (IMF) between 0.1 and 100 
Mq and sub-solar metallicity (0.2 Zq) so that compar- 
isons to existing works are more straightforward. Our 
selection of sub-solar metallicity models is based on the 
observed trends at high redshifts (e.g. Maiolino et al 
2008) as well as on the direct observations of extremely 
blue UV slopes of the most recently found z ^ 7 sources 
(Bouwens et al. 2010b, Labbe et al. 2010). However, 
since other IMFs (or metallicity models) fit our obser- 
vations just as well, we will remark on how the results 
change if we adopt a different IMF or metallicity. For 
simplicity, we assume a constant star formation (CSF) 
rate when modeling the star formation histories of the 
galaxies in our sample. This assumption seems prefer- 
able to an exponentially-decaying SF history (which im- 
plies increasing SFRs with redshift), as there is currently 
no evidence that UV-selected sa mples form stars at a 
faster rate at ear lier epochs fe.g.. [Papovich et al.ll200ll 
IStark et al.ll2009l and others). 

Finally, given the challenges in constraining both age 
and reddening based upon the photometric informa- 
tion available for individual sources, we will assume 
that galaxies in our sample show negligible dust extinc- 
tion when doing the stellar population modeling. We 
have good reasons for making this assumption, z > 
5 galaxies have been found to have very blue J/F- 
continuum slopes and thus little dust extinctio n (e.g., 
Lehnert fc Bremer 2003: Stanwav ct al. 2005; Ya n et al.l 



20051: IBouwe^IsinLl[2006L i2009ai) . As we will show later, 
a similarly steep ?7T^-continuum slope /3 (—2.4 ± 0.4) is 
found for the mean SED of our sample (§6.2) - again sug- 
gesting minimal dust extinction. For individual sources, 
however, it is difficult for us to obtain useful constraints 
on the J/y-continuum slope and hence dust extinction. 
While we have high quality iJieo-band fluxes for our can- 
didates, the other fluxes we have which probe the UV- 
continuum are not adequate: the Jno-band fluxes we 
have available depend significantly on the redshift of the 
source (the Jno-band extends to ^8000 A) and the K- 
band fluxes (available on fewer sources) are more un- 
certain in general (due to the shallower nature of the 
i^-band imaging data). We have tested the impact of al- 
lowing modest amounts of extinction (Ay < 0.5) to the 
models and found that the main results (SMD, SSFR, 
ages) are unchanged or at least consistent within the un- 
certainties. 

Within these general spcciflcations, wc explored a wide 
variety of different parameters (redshift. Age) in model- 
ing the observed photometry of each 2^7 candidate: 



log {Age[Myr]): 



= 4.0 - 11.0 (steps of 0.01) 
= 7.5 - Agcmax (steps of 0.01) 



where Agcmax is the age of the universe at the cor- 
responding redshift so that the models avoid solutions 
in which the populations are older than the universe. 
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Fig. 3. ^ — Observed SEDs and sy nthetic stellar population model fits. The error bars and upper limits are la. The models shown here 
are from lBruzual Hi Chariot II2003I ) with a Salpeter IMF and 0.2 Zq. The red curve represents the best fit models without extinction and 
the blue curve shows the eiTect of imposing a maximal extinction of Ay = 0.5. The open circles shown for the 3.6 iim flux measurements 
of objects GNS-zD2, GNS-zD5, and HDFN-3654-1216 correspond to points with poor photometry. 
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TABLE 3 

Summary of BC03 Models fit parameters. 0.2 Zq, Salpeter IMF, constant SFR, Av = 0. 



ID 


'^phot 


Mass 


Age^'' 




SFR 


SSFR'' 


Luv U-V Hi60 - 3.6 x?ed 






[Myr] 






[Gyr-i] 





UDF-640-1417 

UDF-983-964 

UDF-387-1125 

UDF-3244-4727 

GNS-zDl 

GNS-ZD2 

GNS-ZD3 

GNS-ZD4 

GNS-ZD5 

CDFS-3225-4627 

HDFN-3654-1216 



6 9+°-^ 

"■^-0.1 


c c+0.3 


379 


13.2 




2 0+° =* 
^•'-'-0.1 


5.4 


0.5 


0.9 


1.7 


7 q+0.4 
' ■'^-0.3 


9 9+2.0 
^•^-1.5 


173 


9.0 






3.7 


0.1 


0.6 


0.4 


7 -1+2.0 




19 


4.8 


K 9+5.0 


9cr ^+12.3 


2.0 


-0.7 


-0.8 


0.3 


7 Q+0.8 


2.8tO:5 


315 


6.8 


5 4+2-0 
^•^-1.2 


9 4+6.9 
^■^-0.4 


2.8 


0.4 


0.4 


0.9 


7 9+0.2 
'■^-0.2 




362 


15.8 


12.6+^-1 


2 l+'^-i 


6.5 


0.5 


1.3 


1.5 


7 1+1.5 




251 


7.3 


5 9+'5 * 




3.0 


0.3 


0.3 


0.4 


7 O+0.9 
'■■^-0.4 


4.2«-j 


354 


8.8 


7 -1+2.0 
' -^-1.1 


2 1+^-2 


3.6 


0.4 


1.0 


0.3 


7 9+0.4 
'■^-0.2 




362 


13.9 


11 2+^-'' 


9 1+0.3 

''■^-o.i 


5.7 


0.5 


1.1 


0.8 


7 O+0.2 
' ■■'-0.2 


i2.3j;«? 


354 


26.1 


20.9+J-4 


9 -1+0.4 
^■-^-0.0 


10.7 


0.4 


0.7 


2.3 


7 1+1.5 


3 5+^-'' 


281 


9.2 


7.4l?i 




3.8 


0.3 


0.3 


0.6 


"■■J-0.2 


6^931 


426 


12.4 




-■-■^-o.i 


5.1 


0.5 


0.7 


0.8 


7 o + O.l 

'■•'-o.o 


6.3i;;;i 


354 


13.4 


10.7i«J 


9 -1+0.0 
^•^-0.0 


5.5 


0.4 


0.7 


3.3 



Mean SED'= 



Note. — Best fit parameters, 68% confidence intervals and corresponding for lBruzual fc CharlotI II2003I ) fits with subsolar 
metallicity (0.2Zq) and Salpeter IMF between 0.1 and 100 Mq. We obtain a redshift for the sample of 7.2±0.5. The models 
suggest that in some of these galaxies most of the stars were born at considerably earlier times (z < 10), well into the epoch of 
reionization. The masses of these sources range 0.2-12x10® Af©- We have restricted the dust extinction to zero, consistent with 
the trend observed a t these high redshit s. SFRs from the models are consistent with the ones derived from the extrapolated 
^1500 using the usual lMadau et al.lll99^ formula. These SFRs are somewhat high, ~ lOA/0yr~^ typically. The mean SSFRs of 
the sample is 2.4it0.6 Gyr~^, with an outlier corresponding to the youngest model. The estimated best fit ages do not change 
considerably if we consider either Chabrier IMF or solar metallicity models. As expected for a Chabrier IMF, the masses 
derived are a factor ~0.55 lower and so are the SFRs. Solar metallicity models produce masses ^ 10% larger. The maximal 
extinction model with Ay = 0.5 produces ~ 45% lower ages and ~ 45% higher masses (with the consequent increase in SFR) 
with more scatter. 

Age™ corresponds to the SFH-weighted ages. In the case of CSF models, this simply corresponds to half the time since the 
onset of star formation. The typical uncertainties in this quantity are substantial (1^20 '^y)- 

The SSFR (= SFR/Mass) here is derived from the extrapolated Lisoo luminosity and the masses from the CSF models. 

These are not the mean values of the best fit parameters derived for the sample but rather the parameters derived from the 
modeling of the mean SED described in §6.2 (see figure|4]|. 
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To derive 68%, 95%, and 98% confidence intervals for 
the above parameters, we ran a number of simulations 
where we added photometric scatter (noise) to the ob- 
served fluxes and then used the results to determine the 
threshold values that encompass these confidence in- 
tervals. We performed the above calculations with the 
stellar po pulation modeling code named FAST (see ap- 
pendix in lKriek et al]|2009[) . 

Note that in stellar population modeling, we do not use 
the photometric redshift estimates fi'om the previous sec- 
tion. This was to avoid additional complications in the 
definition of the confidence intervals (as the templates 
are slightly different). However, we found that the two 
redshift estimates are consistent (the typical discrepan- 
cies arc rms Az ^ 0.04). 

In figurelSl the observed SEDs are presented along with 
the best fit models with no extinction. The correspond- 
ing properties we derive for these galaxies using these 
models are detailed in tabled The ages in that table cor- 
respond to the average age of the stellar population, i.e., 

^J^totiil I '^lookback SFR{tiookback)dtiookback (=agCu,). In 

the case of constant SFR models, this is simply equal 
to one half of the time elapsed since the onset of star 
formation. Similar to what others have found at some- 
what lower redshifts, these models indicate the presence 
of quite massive systems very early in the universe. We 
find masses in the range of 0.2 — 12 x lO^M©. These 
models also show quite large ages that place the forma- 
tion of most of the stars up to 380 Myrs earlier with 
typical values in the order of 300 Myr. The uncertainties 
in the derived age, however, are quite substantial, typ- 
ically ^220 (68% confidence intervals). From table 
[3l it can be noticed that the values of age^, imply that 
some galaxies have been forming stars for times compa- 
rable to the age of the Universe at this redshift. We have 
preferred not to impose any arbitrary restrictions on the 
time of onset of star formation because any meaningful 
redshift constraint (e.g. z<100) would not meaningfully 
restrict the ages (16 Myr for z <100), particularly when 
compared to the uncertainties associated with the esti- 
mation, typically ^^20 Myr. Such constraint implies in- 
significant fractions of stellar mass assembled at extreme 
redshifts. 

6.1. Parameter dependencies 

In the above stellar population modeling, we adopted a 
Salpetcr IMF and assumed sub-solar mctallicity (0.2Zo). 
We also explored the effects of using a Chabricr IMF 
and of varying the metallicity of the models. We find 
that while distribution of best fit redshifts and ages is 
unchanged, the derived masses are ^45% smaller if we 
use Chabrier IMF (instead of Salpeter), and 10% larger 
if we consider models with solar metallicity (instead of 
0.2Zq). For both metallicities (0.2 Zq and Zq) and both 
IMFs (Chabrier and Salpeter), we obtain reasonable 
fit results, so there is no reason to prefer one IMF or 
metallicity over the others. 

We have also assumed that our sources suffer from 
minimal dust extinction (both due to the very blue UV- 
continuum slope /3 measured from the mean SED here 
[§6.3] and due t o the blue t/F-continuum slopes ob- 
served at z > 5 : 'Lclmert & Brcmcr'2003'; 'Stanwa v et al.l 
l2005HYan eral..,2005;.Bouwens et al...2006,. 2009a,) . How- 



ever, this extinction is not very well constrained for in- 
dividual z ^ 7 galaxies (which lack strong constraints on 
their C/l^-continuum slopes), and so it is worthwhile to 
mention how larger values of the dust extinction would 
affect the ages and stellar masses derived from our mod- 
eling. Imposing an extinction of Ay = 0.5 (following 
Calzetti et al. 2000 - we consider this a safe upper limit 
based on the previously mentioned studies), yields best 
fit models that are ~ 40% more massive (to compensate 
for the dimming) but ~ 45% younger (which prevents 
their H — [3.6] colors from becoming too red to match 
the observations) with the consequent increase in SFR. 
We also find reasonable for these models. 

Finally, we must remember that we adopted a spe- 
cific form for the star-formation history - supposing the 
SFR for each was constant in time. We could easily 
have adopted other star formation histories (instanta- 
neous burst, exponentially decaying) in fitting the ob- 
served SEDs and found acceptable results. To deter- 
mine the approximate effect of the star formation history 
on our derived parameters, we also considered exponen- 
tially decaying e^*/"^ histories and instantaneous bursts 
in modeling our sources. In general we found larger ages 
and masses for histories with larger t's (where r for a 
constant SFR model is of course 00 and for an instan- 
taneous burst) though all assumed histories produced ac- 
ceptable fits. In that sense, the unrealistic instantaneous 
burst models provide a lower limit for the age^, of the 
sample that we find to be 80 Myr. Using a more reason- 
able approximation like r = 100 Myr this number goes 
up to 170 Myr. This latest result would imply that we 
are observing a quiescent population right after the main 
star formation episode is over, which seems unlikely. 

While we note these dependencies, we will use our fit 
results assuming no dust extinction, a Salpeter IMF, sub- 
solar metallicity (0.2 Zq), and a constant SFR when de- 
riving results in subsequent sections and in particular to 
estimate the Stellar Mass Density (SMD) of the universe. 
Relaxing our assumption that there is no dust extinction 
{Ay=0) to Ay < 0.5 produces no significant change in 
the derived quantities. 

6.2. Average Spectral Energy Distribution 

One significant challenge in modeling the stellar popu- 
lations of z-dropouts in our sample is the faintness of the 
sources and therefore the still sizeable uncertainties on 
the fluxes we derive. Consequently, it becomes difficult 
for us to obtain tight constraints on the model parame- 
ters - like dust or age - for individual sources. 

We can obtain much tighter constraints on the proper- 
ties of our z ~ 7 galaxy candidates by averaging the 
measured fluxes for the sources and deriving a mean 
spectral energy distribution, per unit wavelength. This 
is particularly valuable for a determination of the UV- 
continuum slope /3 for the sample, since this slope is 
constrained from the H and i^T^-band fluxes and since 
the i^-band flux is only poorly constrained for individ- 
ual sources. To derive this mean SED, we first normal- 
ize all the sources to the average iJigo-band flux and 
then take a weighted mean of the fluxes of all sources 
(< X )/(E The mean K-hand flux 

is derived from the 6 sources where we have deep K- 
band data. The sources with poor photometry (GNS- 
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estimate of the stellar mass density in §7. 
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Fig. 4. — The average observed SED (black circles) for galax- 
ies derived from the 11 sources in our z-dropout selection. Error 
bars and upper limits are at la. The fluxes of each source in our 
z ~ 7 2-dropout selection was scaled so that their iTigo-band flux 
matched the average ffigo-band flux for the sample, after which 
the rescaled fluxes from the entire sample were averaged. Flux 
measurements which were poor (due to difficulties in precisely sub- 
tracting a nearby neighbor, as for the 3.6fiin flux measurements of 
GNS-ZD2, GNS-zDS, and HDFN-3654-1216) were not included in 
the average. Note, the average A'-band flux we derive here was 
only determined from the 6 sources for which we had deep X-band 
data. Stellar population models were fit to this average observed 
SED using the same technique that what was used for individual 
sources (model shown in solid line). The parameters of the fit are 
presented in table|3] The average zgso-band flux (open circle) was 
not included in the flt - since it would not make sense to include 
both the 2850 and Jno-band flux (both of which have a different 
dependence on redshift) in the fits. This average SED shows no 
detection in the optical B, V, or i bands and a very large break 
(>3 mag) between the optical and J bands - suggesting that the 
majority of sources in our sample do in fact correspond to z ~ 7 
galaxies. The extremely blue H — K color, which traces the UV 
continuum, imply that the dust reddening must be very low — con- 
sistent with our assumptions in modeling individual z-dropouts 
(§6). The pronounced break (~1 mag) between the K and the 3.6 
fim bands suggests the presence of a IJalmor break, indicating that 
the typical z ~ 7 galaxy has experienced several previous genera- 
tions of star formation. We have used the derived M/L ratio from 
this mean SED to make a simplified estimate of the stellar mass 
density of the universe at z ~ 7 of 4.5x10^ AfQAfpc-3 (see §7). 

zD2, GNS-zDS and HDFN-3654-1216 in 3.6fim) were not 
considered when taking the mean. The mean SED is 
presented in Tableland Figure [3] The C/F-continuum 
slope /3 we estimate from the H and K photometry is 
-2.4 ±0.4. 

We perform the stellar population modeling for this 
source in the same way as for the other sources and in- 
clude the estimated properties in Table [H When mod- 
eling, however, the Zs5o flux has not been included in 
the fits to minimize the influence of our lowest redshift 
sources. The best-fit model is also compared with the 
mean SED in Figure |3l If we allow the dust reddening to 
be non-zero and include that in the fits, we find a best 
fit Ay = 0.4. This provides support for our assumption 
in §6 that our z ~ 7 galaxy candidates are largely dust 
free. Using the mass of the best fit to the mean SED 
(6.3x10^ Mq) and the extrapolated luminosity at 1500 
A (5.5x 10^° Lq) we infer a mass to light ratio of Ad/Luv 
= 0.12 Mq/Lq, and also a M/Ly = 0.28 Mq/Lq (from 
Ly = 2.2 X 10^° M0 measured from the best fit model 
at 5500 A). We will use these M/L ratios to provide one 
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Fig. 5. — The SSFR measured from our data at z ~ 7 com- 
pared to the values derived from the data presented by others at 
a constant stellar mass of 5xlO^M0 (corresponding to the me- 
dian of the pres ent sample). Our es timates of the SSFR are based 
Noeske et aI7 f2007). D'adcfi'crar. (2007') (in good 



on results from ^^^^^ 

agreement with lPapovich et al . 2001). and Stark et al. (200jJ). We 
estimate that the typical errors at z< 7 (dashed lines) are ~ 0.3 
dex. The SSFR seems to be remarkably constant at 2 Gyr"'^ be- 
tween z ~ 2 — 7 suggesting that the star formation - mass relation 
does not evolve strongly laetween z ~ 2 — 7. The drop observed 
at z < 2, however, indicates that some physical process might be 
inhibiting star formation. This result suggests that star formation 
in galaxies at z > 2 follows somewhat different principles than for 
galaxies at z < 2. 



6.3. Specific Star Formation Rate 

A key quantity in considering the build-up of stars 
within a galaxy is the specific star formation rate - sim- 
ilar to the b parameter which was more frequently used 
in the past to characterize a galaxy's star-formation his- 
tory. It is the star formation rate within a galaxy divided 
by its stellar mass - or equivalently the fraction of the 
stellar mass in a galaxy that forms per unit time. As 
such, the SSFR provides us with a useful way of think- 
ing about galaxy growth over cosmic time - and so it is 
not surprising that it has been estimated out to z ~ 6. 

Here we make use of our best fit synthetic models to es- 
timate Li5oo and then use the usual formula from Madau 
et al. (1998) to estimate the SFR. We combine this with 
the masses of the best fit models to obtain the specific 
SFR (see caption in table 3). We find values of 1.8 - 
4.1 Gyr~^ for this quantity across our sample, with a 
median value of 2.4 ± 0.6 Gyr~^ (the outlier at ^--^20 
Gyr"-'^ is undetected in the mid-IR so the constraints 
are poor). To put this value in context, it makes sense 
for us to compare our derived SSFRs with the values 
at lower redshift. Given that the SSFRs can depend 
somewhat on stellar mass, we compare our results to the 
same median mass as we fin d in our sample, 5 x 10^ Mq. 
From the data presented in lStark et al.l (|2009f ). we find 



their 



4 , 5, and 6 sam- 



2.1, 2.0, a nd 2.0 Gyr" . 

pies ; from [Papovich etHI (j2001[) . ISawicki et all (|2007f ) 
and lDaddi et al.l (|2007l) w e find - 2 Gyr"^ at z - 2; and 
from iNoeske et all ()2007D we estim ate 0.3-1.2 Gyr~^ 
for z ~ 0.2 - 1 samples. At z ~ 2, iReddv et all (|2006( ) 
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find a much higher value of ~ 10 Gyr~^ but at that mass 
it is only based on MIPS detected sources. Obviously, 
any MlPS-detected sample would be biased to include 
only those sources with substantial enough SFR to show 
MIPS detections, and hence probably is not representa- 
tive for the lO^'^Af0 population. The results are shown 
in Figure [5] 

From z ^ 7 to z ^ 2, there is little apparent evolution 
in the SSFR and the present re sults provide a con tinu- 
ation of the tren d delineated by 'St ark et al.l (I2007D and 
lYan et all (|2006D from the Jleddy et al.l (|2006ir points. 
However, from z ^ 2 to z ^ the SSFR shows a rapid 
decrease. This suggests that SFR at z > 2 mostly pro- 
ceeds in a largely similar way, but that at z < 2 there 
must be some physical processes inhibiting SFR. Simi- 
lar to many other c omparisons of merit, e.g., evolution 
of M* with redshift (iBouwens et alll2006l: lYoshida et al.l 
120061: IBouwens et al.l l2007D . we find that SF in galaxies 
at z > 2 seems to follow somewhat different principles 
than for galaxies at z < 2. 

7. STELLAR MASS DENSITY 

One of the most fundamental quantities we can try to 
infer from z ^ 7 selections and the present stellar pop- 
ulation modeling is the stellar mass density. The stellar 
mass density tells us how much star formation occurred 
in the universe to the point of observation, and there- 
fore provides us with a very powerful constraint on early 
galaxy formation. 

7.1. Selection Volumes 

An essential step in determining the stellar mass den- 
sity from our z ~ 7 z-dropout selection is obtaining an 
accurate estimate of the selection volume. This requires 
that we model the selection of dropout galaxies from all 
six of the Bouwens et al. (2009, in preparation) search 
fields and estimate the effective volume we are able to 
search versus i/igo-band magnitude. The selection vol- 
umes are calculated by adding artificial sources to our 
search fields and then attempting to reselect them as 
z ~ 7 z-dropouts according to the criterion described in 
§2. The artificial sources are assumed to have a mean 
J/y-continuum slope /3 of —2, consi stent with the ob- 
served trends a 4 < z < 6 (jBouwens e t al. 20093) ■ Their 
pixel-by-pixel surface brightness profiles are identical to 
those of random z ^ 4 i3-dropouts from the HUDF 
(jBouwens et al.|[2007|) of similar brightness but their sizes 
have been rescaled as (1 -I- z)~^, following the observed 
size trends with redshift at z > 2 fo csch ct al.. 2009d|; 
iFerguson etldl 120041: [Bouwens et alTl2004bl) . We have 
not included the possible contribution from Lyman al- 
pha emission to the SEDs. A more detailed description 
of these simulations can be f ound in Bouwens et ah (2009, 
in preparation: but see also IBouwens et al]l2008[ ). 

We estimate the following total search volumes as a 
function of -ffigo-band magnitude for the entire Bouwens 
et al. (2009, in preparation ) z ^ 7 z-dropout selection 
(from both the HST NICMOS and ground-based near-IR 
data): 



VeffiH. 



160) 



33 X 10" Mpc^ if 25.3 < Hiea < 25.8 

16 X W Mpc^ if 25.8 < Hiea < 26.3 

6.0 X 10''Mpc^ if 26.3 < Hiea < 27.1 

1.2 X 10-* Mpc^ if 27.1 < Hieo < 27.5 



7.2. Stellar Mass Density Determinations 

We proceed now to estimate the stellar mass density 
of the sources we find in these search volumes. We do so 
in the three different ways we detail below: 

Direct Approach: Here we estimate the stellar mass 
density by simply suming over the expected mass den- 
sity to come from each source. The expected mass den- 
sity is simply the estimated stellar mass for a source 
multiplied by the the likelihood it is not a contaminant 
(90% for sources in the HUDF and 80% otherwise) di- 
vided by the selection volume above. The stellar mass 
density we derive by summing over the 11 sources is 
5.7x10^ Mq Mpc~^. We estimate the uncertainties in 
the mass density by bootstrap resampling (detailed be- 
low). This approach has the advantage of being very 
direct and even-handedly including all candidates in our 
selection in the estimate. The disadvantage, of course, is 
that this estimate may be affected by the mass estimates 
of the three z-dropouts with perhaps unreliable 3.6 fim 
fluxes. 

Mean SED Approach: For our second estimate, we use 
our average SED from §6.2 to derive the mean M/Luv 
ratio for sources in our sample - reasoning that this M/ L 
ratio is much more accurately known than any individ- 
ual Al/L ratio in our sample. For our second estimate, 
we use our average SED from §6.2 to derive the mean 
M / Ljjv ratio for sources in our sample - reasoning that 
this M / L ratio is much more accurately known than any 
individual M/ L ratio in our sample. If the photometric 
redshifts are accurate (z ~ 7.2 ± 0.5), then the limiting 
depth of our search co rresponds to ~ Mjjvab = —20. 
We have then used the IBouwens et all (|2008f ) UV LF at 
z ^ 7 to estimate the UV Luminosity density integrated 
to that depth and multiplied it by the MjLuy ratio 
to obtain a stellar mass density of 4.5x10^ Mq Mpc""^. 
Again for the uncertainties, we rely on a bootstrap re- 
sampling procedure (below). 

Random M/L Approach: For our final estimate of 
the stellar mass density, we ran a monte-carlo simula- 
tion where we match up the 11 galaxies in our z 7 
z-dropout sample with randomly sampled M/L ratios 
from the 8 galaxies with reliable IRAC fluxes. We then 
divide the masses by the search volumes that corre- 
spond to their UV luminosities. After repeating the 
match up process 100000 times we find a median value of 
6.6x10^ Mq Mpc~'^. The estimate of the uncertainties 
based on bootstrap resampling is described below. 

Uncertainty Estimates: To estimate the uncertainties 
in our estimate of the stellar mass density at z ^ 7, we 
must fold in the many uncertainties that contribute to 
this density, including uncertainties in the sampling of 
the LF, the M/L ratios of the galaxies we uncover in our 
search, whether any individual source is a contaminant, 
and finally the mass estimates themselves of candidates 
in our sample. 

The simplest way of including all these uncertainties 
in our final estimate is to run a Monte-Carlo simulation. 
For each simulation, we iterate over all 11 candidates in 
our sample and include Pi (x) number of sources in that 
trial with a UV luminosity equal to that candidate - 
where Pi {x) is a Poisson distribution with mean equal to 
1. Then, we run over all the sources and give each source 
10-20% chance of being thrown out (to account for uncer- 



13 



10° 



^10^ 

■4-J 

0) 
Q 



4-' 



10 



• Mi5„„<-20, Stark et al. 09 

• Labbeetal. 06 
Current sample 



+ 



3.5 4.0 4.5 5.0 



5.5 6.0 
Redshift 



6.5 7.0 7.5 8.0 



Fig. 6. — The stellar mass density as a function of redshift. Our 
estimate of the stellar mass density at z ~ 7 is shown with the red 
star and considers the contributions from galaxies with absolute 
magnitudes Muy^ab < —20. The horizontal error bars show the 
approximate width of our samples in redshift s pace. Also shown 
are the stellar mass density determinations of IStark et al.l 120091 
at 2 ~ 4, 2 ~ 5, and z ~ 6 (b lue points). The green point at 
2 ~ 7 is from lLabbe et ani2006l . That estimate w as derived from 
two 2 ~ 7 candidates in the HUDF identified by IBouwens et al.l 
l|2004d ). Those candidates are also included in this study. The 
limiting UV luminosity probed by our sample is comparable to the 
points at lower redshift. 

tainties in the contamination fraction). Next, we assign 
a mass to each of the objects in the particular reahzation 
by drawing a random M/L^y from the observed values. 
When doing this, we also include individual uncertain- 
ties in the M/L^/y determinations, specifically, we make 
a weighted choice of a mass from the distribution asso- 
ciated to a particular source (with the weights derived 
from the as was already described for the determina- 
tion of the ages). Finally, we divide each source by the 
selection volume appropriate to its UV luminosity and 
sum the sources to calculate the stellar mass density for 
a given trial. We repeated the simulation 10000 times to 
ensure that our results were not limited by the number 
of trials. After sorting the distribution, we found that 
the 68% upper and lower limits on the stellar mass den- 
sity were 3.3 x lO'^ Mq Mpc~^ and 1.2 x lO*' Mq Mpc"^, 
respectively. 

Summary: Above we derive three different estimates 
of the stellar mass densities with uncertainties. The es- 
timates are 5.7, 4.5, and 6.613:3 x lO'^ Mq Mpc"^. Ah 
these estimates are consistent with each other but we pre- 
fer the random M/L approach because this one should be 
less affected by possible systematic errors in the mass de- 
rived from poor photometry. We present this stellar mass 
density in Figure [6] and show the previou s determina- 
tions o f ISterk et alT(|2009[ ) at z 4 - 6 and lLabbe et all 
()2006D at z ^ 7 for context. We discuss differences be- 
tween these stellar mass density determintions and the 
observed trends in §8 and §9. We should also note here 
that allowing moderate extinction {Ay < 0.5) produces 
SMD measurements that are fully consistent with the 
previous values (within the uncertainties). 

7.3. SFR Density Determinations 
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Fig. 7. — The approximate star formation rate density inferred at 
2 ~ 8.5 (solid red circle) by combining the age constraints we have 
on 2 ~ 7 galaxies with the stellar mass density we estimate (§7.3). 
For our estimated SFR density, we use the stellar mass density we 
derive at 2 ~ 7 divided by the average age of our sample. Included 
on this figure are the dust-corrected S FR density determin ations 
at 2 ~ 7 from 2 ~ 7 2-dropout search I IBouwens et ai.ll2008l : open 
red circle), the IBouwens et al.l 1120 0"'^) de termination at 2 ~ 4 — 6 
(open red circles), the IReddvfc^ teidel (2009) determinations at 
2 ~2-3 (green crosses), and the [Schiminovich et al. (200.5) deter- 
minations at z < 2 (b lack hexagons). Dust corrections are as in 
IBouwens eFall f2009ai V 

The advantage of the current stellar population mod- 
eling is that it permits to estimate the SFR in our can- 
didates at even earlier times. Combining our age con- 
straints with the estimated SMD we can place limits on 
the SFR density at even earlier times. 

In terms of the previous history, the constant SF mod- 
els imply ages of the galaxies consistent with them being 
in place at around z ^ 10, which combined with their 
assembled masses imply a very simple estimate of the 
average SFRD between 7 < 2; < 10. The amount of time 
elapsed corresponds to 300 Myr and the stellar mass as- 
sembled in that amount of time is 3.3 x IO^MqMpc"'^. 
This implies a SFRD of 1.1 x 10"^ Mq yr"! Mpc"^ (see 
Figure [7]) . An extreme approach that would maximize 
the SFRD comes from the single burst models. Simply 
dividing the total masses (6.6 x IO^MqMpc"'^) by the 
ages of the sources (80 Myr for the single burst models) 
yields an average SFRD of 8 x 10"^ Mq yr'^ Mpc"^ in 
the previous 80 Myr. 

8. RELIABILITY OF CURRENT RESULTS 

8.1. Comparison with Previous Photometry 

In general, our optical to near-IR photometry is con- 
sistent with the photometry presented in Bouwens et 
al. (2009, in preparation) although a systematic off- 
set of ^ 0.2 magnitudes is present due to the fact that 
the fluxes presented there were measured in a somewhat 
smaller aperture and no aperture corrections were ap- 
plied. We also find excellent agreement between the 
mid-IR IRAC fluxes me asured for our HUDF z-dropouts 
and those presented in iLabbe et al.l (|2006f l. Although 
this might not be surprising due to our use of the same 
technique for doing photometry, our modeling of the flux 
from neighboring sources (and subsequent removal of this 
flux) is completely independent. This illustrates the ro- 
bustness of our method for doing photometry. An in- 
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dependent test of the quality of photometry can be ob- 
tained by comparing flux measurements for sources with 
deep IRAC observations taken in both GOODS epochs 
(rotated by 180 degrees). In general, we observe excellent 
agreement between the two measurements for the five 
sources with two epoch data (the four UDF z-dropouts 
and CDFS-3225-4627) - suggesting that systematics are 
minimal. The only exception to this is for the 3.6 fj,m 
measurement for GNS-zD2, GNS-zD5, and HDFN-3654- 
1216 where there are bright nearby neighbors. 

8.2. Comparison with Previous Estimates of the Stellar 

Mass Density at z ^ 7 

Previously, iLabbe et al.l (j2006() made an estimate of 
the stellar mass density at z ~ 7 based upon a small 
(4 galaxy) HUDF z-dropout sample. They estimated 
a stellar mass density of 1.6lj j x 10*^ Mq Mpc^^ to 
>0.3 L*^^. Since we adopt a similar limiting luminosity, 
we can make a direct comparison with the stellar mass 
density estimated here. We find a fiducial value that is 
a about half the one estimated by ILabbe et al.l ()2006D . 
This is mostly due to our different choice of star for- 
mation histories (we favor constant SF histories versus 
the average between SSP and constant SFR used in that 
work). These estimates are fully consistent within the 
uncertainties. 

This quantity was also estimated bv lStark etldl (|2009f ) 
at z ~ 4 — 6 for sources in the GOODS fields to similar 
depth (see figure O. A simple calculation shows that the 
observed growth in mass between z ~ 7 and z 6 is 
consistent wit h the observed SFRD derived from the UV 
LF studies of iBouwens et al.l (l2008f ). 

8.3. How Significant Is Crowding for Current Samples 

Because of the broad PSF in IRAC data, crowding is 
considered to be a potentially significant concern in do- 
ing photometry on faint sources, particularly when these 
sources are nearby bright foreground galaxies. In fact, in 
many studies, it is thought that perhaps ^--^50% of faint 
sources are sufficiently close to their neighbors that IRAC 
photometry is impossible. What do we find here? 

We attempted to do IRAC photometry on all 11 z- 
dropouts in the Bouwcns et al. (2009, in preparation), 
without excluding any sources due to crowding issues. 
Of the 11 sources presented here, 5 suffered significant 
blending with nearby neigh bors. However, as a result of 
our deblending technique (jLabbe et al.l l200"6l ). we were 
able to recover reliable fluxes for all 11 in the 4.5 fj,m 
band and 8 of 11 (73%) in the 3.6 /im band (exclud- 
ing GNS-ZD2, GNS-zDS, and HDFN-3654-1216 which 
showed strong residuals from the neighbors after the sub- 
traction process - implying large systematic errors in 
those two cases). This agrees with Monte-Carlo experi- 
ments that we performed that suggest that photometry 
is possible for >80% of faint sources and that the largest 
errors in the recovered fluxes should be roughly a factor 
of 2. 

9. DISCUSSION 

9.1. Stellar mass growth during the first Gyr of the 
Universe. 

Despite the large uncertainties in the derived indi- 
vidual masses, it seems clear now that quite massive 



(> 10^° -^^o) systems with evolved stellar populations 
were already present in the universe at very early times 
(z ^ 7 or at least z ^ 6). These very massive systems 
likely correspond with the most massive Dark Matter 
Haloes (DMH) and arc predicted to exist but in low num- 
bers in the standard Press Schechtcr formalism. It is 
hard to run simulations that can probe this very mas- 
sive end of the Mass function since to obtain them in 
larger numbers it would be necessary to run simulations 
with quite large volumes. iDave etldl (|2006| ) ran an SPH 
simulation with a comoving volume of 10^ Mpc^ which 
they find to be adequate to probe the stellar mass range 
10^'^ — 10^° Mq. At the median mass of our sample of 
5x 10^ Mq, they find a number density of ^ lO''^ Mpc"^. 
More recently, Choi & Nagamine (2009, in preparation), 
ran a simulation in a larger box of ^ 3 x 10^ Mpc"^ and 
find the number density at 5 x 10^ M(t, to be a pprox - 
imately one half that predicted by IDave et~dl\ ()2006[) . 
Given the uncertainties, these predictions agree quite 
well with the density of such objects that we find of 
1.6 X 10~^ Mpc~^. Although our sample may still suf- 
fer from small number statistics, the fact is that such 
massive objects are found to be fairly common at this 
early epochs and their moderately evolved SEDs show 
that their stellar populations arc not pristine but were 
partly formed at higher redshifts. 

9.2. Star Formation Histories of high redshift galaxies. 

It is well known that the results (e.g., derived masses, 
SFRs, etc.) of stehar population analyses can depend 
significantly on the functional form one assumes for the 
star formation history. Some care, therefore, needs to be 
given to the parameterization of these histories to ensure 
that the conclusions drawn do not depend too much on 
artifacts of this parameterization. 

To illustrate, there are 3 different model parameteriza- 
tions of the star formation history in common use in the 
literature: instantaneous burst models, exponentially de- 
cay star formation models, and constant star formation 
models. Instantaneous burst models give younger ages 
and lower masses than exponential decay models which 
in turn give younger ages and lower masses than constant 
star formation models. The instantaneous burst models 
do not seem realistic but are useful to set lower limits on 
the derived ages and masses. Exponential dec ay models 
are the most popular i n the literature ( e.g. lYan et aP 
I200I lEvles et alJ [20071 : iStark et al.l [20091 1 - perhaps be- 
cause of their versatility in modeling a wide range in star 
formation histories - but give rise to a troubling predic- 
tion, namely, that the star formation rates of galaxies are 
larger in the past. This prediction is troubling because it 
contradicts both the observed and predicted trend that 
the SFR density of the universe increases with cosmic 
time. 

Given these concerns, we prefer to model the stellar 
population of high-redshift galaxies with constant star 
formation histories. The reason for this preference is 
as follows: first, constant star formation models do not 
naturally predict that the SFR density will be greater 
at early times, as exponentially decaying models do. 
Second, constant star formation models do not predict 
a large population of UV luminous sources at earlier 
times as found from exponentially-decaying models (e.g.. 
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TABLE 4 

Key results derived from z ^7 sample. 



Quantity 



Value 



Rcdshifts 
Masses 
M/Ljjy ratio 
Minimum Age'' 
Average Age 
(7y-continuum Slope f} 

Specific SFR 
Mass Density (Direct) 
Mass Density (M/L of Mean SED) 

Mass Density (Random M/L)'' 
Predicted SFR density (at 2 = 9) 



7.2 ±0.5 

0.1-12x10'' Mq 
0.01-0.1 Mq/Lq 
80 Myr 
300 Myr 
-2.4±0.4 
2.4±0.6 Gyr^^ 
5.7x10^ Mq Mpc-3 
4.5x10^ Mq Mpc-3 
tg j X 105 Mq Mpc-3 
0.0011 Mq yr-lMpc-3 



6.6 



^ Prom single burst models. 
^ Our best estimate. 



lYan et al.l[2006l) . This is more consistent with the obser- 
vation that such sources are not found in large numbers 
at z ^ 7. Third, constant star formation models are 
more consistent with the low evolution seen in the spe- 
cific star for mation rate (e.g., Figure [5l here, §6.2, and as 
discussed bv lStark et al.l[2009| ). We realize that some lu- 
minous galaxies may have mass-to- light ratios suggesting 
their SFRs were higher in the past, but we suspect this 
may be a duty cycle issue that is due to feedback, etc. 
These galaxies are simply experiencing a period where 
their SFRs arc less than their norm. 

9.3. Reionization 

As outlined in §7.3, the stellar population modeling 
we do of z '--^ 7 galaxy candidates allow us to estimate 
the SFR density at even earlier times. Having an esti- 
mate of this SFR density is valuable since it allows us 
to assess how much ionizing radiation the star-forming 
population at z > 7 might likely produce ~ and hence 
repose the question about whether z>6 galaxies are ca- 
pable of keeping the universe reionized. From absorp- 
tion studies to bright z ~ 6 QSOs, the process of reion- 
izing t he hydrog en of the universe was just ending at 
z ~ 6 (|Fan et aLlbool iBecker et aLll2001| ) while the five 
year W MAP results suggest it began at least as early as 
z ^11 (iKomatsu et al.ll20"09l) . 

iMadau et al.l ( T999f) presented a prescription to esti- 
mate the critical density of UV radiation necessary to 
reionizc the universe at a given redshift. Upd ated to a 
more current cosmology (jKomatsu et al.ir2009f ). that for- 
mula becomes: 



crit / \ 

Psfr{z) 



0.04 

fesc 



0.0463 



At 



7, for an escape fraction /e. 



Mq yr Mpc 
(1) 

. „„c = 0.1 (e.g. 

iShaplev et al.ll2006f). and f or a clumping factor C = 30 
(but see lPawlik et al.ll2009l who suggest C ~ 6) the value 
of of p'sFR = 0-4 Mq yr~^ Mpc~^. The sources capable 
of producing such radiation, however, remain unknown 
but young O and B stars in early galaxies stand out as 
the most likely candidates given the observed decrease in 
the number density of quasars at high redshifts. 

Based upon the stellar mass density and the mean ages 
we derive in §6 and §7 for the sample, we estimated an 
average SFR density of 0.0011 AIq yr~^ Mpc""^ between 



7 < z < 10 (§7.3). This is more than 2 orders of mag- 
nitude below the SFR density required to reionize the 
universe at z = 7. We can also obtain an upper limit 
to the SFR density by considering the minimal ages ob- 
tained from the single burst models. The mean age of 
80 Myr obtained from these models implies a SFR den- 
sity of 0.008 Mq yr-i Mpc^^ at z = 8 (§7.3). Even for 
this larger value for the SFR density, we are still a fac- 
tor ~50 below that required to reionize the universe. Of 
course, these values are based only on the brightest ob- 
servable sources and so including the mass from galaxies 
of even lower luminosities than our selection limit would 
increase these numbers by a factor of > 2 — 3. In any 
case, it bears mention that we cannot include the contri- 
bution from a population of dust obscured sources that 
we would miss in LBG selections. However, as we have 
already noted in §6.1 and §6.2 there is strong evidence 
that the contribution frora this population is not large 
(e.g., iBouwens et"aIll2009aD . 

10. SUMMARY 

We use the very deep optical, near-IR, and IRAC data 
over and around the two GOODS fields to study the 
properties of a large sample of ^11 z-dropout galaxies at 
z ~ 7. Considered are the ages, stellar masses, redshifts, 
and dust properties of z ^ 7 galaxies. The z-dropout 
candidates were drawn from the very large selection of 
such galaxies from ^80 arcmin^ of deep NICMOS data by 
Bouwens et al. (2009, in preparation). Essential to this 
analysis is the availability of deep Spitzer IRAC mid- 
IR data that give us deep coverage of these sources at 
rest-frame optical wavelengths and hence permits us to 
estimate the stellar mass. T he stellar population mod - 
eling was performed with the iBruzual fc CharlotI ()2003[ ) 
spectral synthesis modeling code with a Salpeter IMF 
and assuming solar metallicity. We have quantified how 
our best-fit properties would change, as a function of the 
assumed IMF and metallicity. 

Our conclusions are as follows: 

• Photometric redshifts place the candidates at 6.2 < 
z < 8.0 with the mean redshift of the sample at 7.2. 
The uncertainties in the redshift for individual can- 
didates are typically Az 0.5 and have been taken 
into account in the derivation of the confidence in- 
tervals for the rest of the properties (see §5 and 
figure E]). 



• We use the results of of IBouwens et al.l (|2009al ) re- 
garding the low extinctions expected at high red- 
shifts to help us better constrain the individual ages 
of the sources. Our CSF model fits yield SFRs that 
are consistent with simple conversion based on the 
Li5oo luminosities. The best-fit ages allow enough 
time to assemble their masses by the redshift of 
observation. This does not seem to be t he case for 
themost massive sources in the work by lYan et aP 
(l2006l) . We argue that this is not likely an effect of 
large extinctions but rather that these sources arc 
the result of some other more complicated history 
probably involving mergers (§9.2). 

• The star formation history weighted ages of the ob- 
served stellar population we derive are in the range 
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of 170 - 420 Myr with a mean of 300 Myr (there is 
an outher with ^ 20 Myr). These ages are consis- 
tent with previous work that place the formation 
of the most massive galaxies at very early epochs. 
In particular, the bulk of the stars in some of these 
sources seem to have formed as early as z ~ 10 (see 
§6 and table [3]). 

• The stellar masses we estimate for individual z ^ 1 
z-dropouts in our sample range from 0.2x10^ Mq 
to 12 xlO^ M0, with a mean for the sample of 
5.1x10^ Mq. The masses we estimate are much 
more well constrained than other quantities - like 
the age - but are nevertheless still uncertain at the 
factor of 2 level (§6 and table [3]). 

• We find that the specific SFRs (SFR/Mass) of 
the sources in the sample range from 1.8 Gyr^^ 
to 25 Gyr~^, with a biweight mean value of 
2.4±0.6 Gyr^-'^. We observe that at comparable 
masses, the specific SFR is surprisingly close to all 
the valu es in the literatur e between z ^ 2 and 

6 (IReddv et al.l [2006t IPapovich et all [2OOII: 



IDaddi et al.l 120071: IStark et all booC The con- 
stancy of this quantity between z ~ 2 — 7, in con- 
trast with its fast decline at z < 2 (|Noeske et al.l 
I2007D . suggests that star formation proceeds in dif- 
ferent ways in these two regimes (see §6.3 and figure 
[5]). 

• Utilizing the estimated selection volumes for the 
Bouwens et al. (2009, in prep) z-dropout search, 
we derive the stellar mass density at z = 7 us- 
ing an approach that randomly samples the M/L 
of the galaxies with the most reliable IRAC pho- 
tometry. We find 6.6;^;3 x 10^ Mq Mpc'^. The 
random M/L approach is preferred as it is less af- 
fected by possible systematic errors in photome- 
try. We tested other approaches, including aver- 
aging the direct fits, yielding very similar answers. 
Our estimate of the global stellar mass assembled is 
consistent with the growth expected based on the 
SFRD measured between 6 < z < 7 and the SMD 
measured to similar depths at z 6 (§7.2). 

• Combining the estimated ages and total assembled 
mass we can derive an average SFR between z ~ 10 
andz-7ofl.lxlO-3Moyr~i Mpc'^. SSP mod- 
els provide us with minimum age estimates that 



in combination with the masses allow us to place 
an upper limit to the SFRD at z '--^ 8. This esti- 
mate is still a factor 50 below the necessary value 
to reionize the univers e at this redshift (following 
the lMadau et al.|[l999l prescription with fesc = 0.1 
and C = 30, eq. 1), in agreement with previous 
works. We emphasize that this estimate is only 
based on the most luminous sources and is proba- 
bly missing most of the UV light, which is being 
produced by sources in the faint end of the LF 
(e.g., B ouwens e t al|[2007t IReddv fc Steidel [200l 
lYan fc Windhorstli2004D . 

The Bouwens et al. (2009, in prep) search based 
on a large area of high quality optical to NIR data 
has provided us with the first sizable sample of can- 
didate sources at z ~7, a step of 200 Myr with re- 
spect to the previo us effor ts (but see also new work by 
Oesch et al.ll2009bl: [B ouwe ns et al.ll2009bl: iMcLure et al.l 
2009bl : IBunker et al.l 120091 : iOuchi et al.ll2009n . Comple- 
menting these data with very deep mid-IR Spitzer IRAC 
imaging we have been able to fit BC03 SSP models to 
estimate the masses and ages of galaxies at z ^ 7. Our 
results suggest that these galaxies had been forming stars 
for >200 Myr and as soon as z~9-10, well into the reion- 
ization epoch. We expect to substantially improve upon 
these results taking advantage of the deep near-IR data 
soon to become available over the HUDF and CDF-South 
GOODS field as a result of the new WFC3 instrument on 
HST. Not only will we substantially increase the number 
of z ~ 7 — 8 galaxies known, but the deeper data and 
improved set of near-IR filters (Y, J, H) will enable us to 
perform much more accurate stellar population modeling 
on individual sources. 
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APPENDIX 



Table [5J in the main text, summarizes the photometry for the sample in the usual AB magnitude system. Since this 
is a z-dropout sample, most of the optical measurements correspond to upper limits which makes it hard to reproduce 
the stellar population modeling results presented in this work. In table [5] of this appendix we provide an equivalent 
table with the actual flux measurements in physical units. 
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TABLE 5 
Summary of Photometry. 



ID 


B435 


""606 


^775 


Z850 


-'no 


^^160 


Ks 


3.6fim 


A.5fim 



UDF-640-1417 

UDF-983-964 

UDF-387-1125 

UDF-3244-4727 

GNS-zDl 

GNS-ZD2 

GNS-zDS 

GNS-ZD4 

GNS-ZD5 

CDFS-3225-4627 

HDFN-3654-1216 



0.18±0.41 
-0.37±0.41 
-0.07±0.37 
-0.41±0.37 
-0.78±1.33 
-0.41±1.04 
-0.55±0.85 
-0.15±1.48 
-0.33±1.07 
0.78±1.26 
-1.26±0.74 



-0.22±0.26 
-0.07±0.26 
0.04±0.22 
-0.22±0.26 
-0.74±0.89 
0.89±0.81 
0.33±0.67 
-0.11±1.07 
O.llitO.81 
0.59±0.92 
0.59±0.67 



-0.07±0.41 
-0.30±0.44 
-0.04±0.37 
-0.74±0.37 
-0.30±1.78 
-2.29±1.59 
-0.33±0.96 
-1.04±1.59 
-1.89±1.22 
-2.77±1.33 
1.11±1.33 



2.88±0.63 
0.04±0.74 
0.85±0.63 
-0.59±0.67 
1.33±2.48 
1.29±1.59 
0.44±1.26 
1.41±1.63 
2.15±1.41 
1.59±1.63 
9.21±1.89 



8.91±1.37 
7.77±1.37 
3.11±1.04 
3.48±1.07 
9.32±2.29 
5.66±1.96 
4.25±2.66 
7.58±2.48 
11.76±3.37 
6.77±2.37 
12.39±2.55 



13.13±2.11 
7.06±1.66 
6.32±1.63 
6.32±1.78 
12.50±2.55 
7.73±2.48 
7.73±2.33 
11.13±3.00 
28.18±3.48 
9.76±2.70 
15.94±2.22 



17.49±3.14 
8.25±4.66 
2.03±4.07 
-1.26±5.51 



24.59±4.07 



10.02±3.40 



30.88±3.51 
12.06±3.51 
2.96±3.55 
9.47±3.55 
40.20±5.07 
9.84±5.33 
18.97±5.07 
31.36±4.77 
52.74±3.66 
13.24±3.55 
30.18±6.99 



16.75±6.88 
8.21±6.88 
0.07±6.99 
7.40±6.92 
26.11±9.69 
11.87±9.80 
15.87±9.65 
26.70±10.21 
25.78±7.06 
14.68±7.06 
19.45±11.80 



Note. — Photometry for the sample. Fluxes arc in units of 10 njy. This table is equivalent to table [2] but here we present the aetual flux 
measurements in the optieal instead of tabulating the la upper limits. 
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